Modified spin-wave study of random antiferromagnetic-ferromagnetic spin chains 
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We study the thermodynamics of one-dimensional quantum spin- 1/2 Heisenberg ferromagnetic 
| system with random antiferromagnetic impurity bonds. In the dilute impurity limit, we generalize 

the modified spin-wave theory for random spin chains, where local chemical potentials for spin-waves 
in ferromagnetic spin segments are introduced to ensure zero magnetization at finite temperature. 
' This approach successfully describes the crossover from behavior of pure one-dimensional ferromag- 

net at high temperatures to a distinct Curie behavior due to randomness at low temperatures. We 
discuss the effects of impurity bond strength and concentration on the crossover and low temperature 
behavior. 
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I. INTRODUCTION 



The renormalization group (RG) approach of WilsonB provided a framework in which to understand the behavior 
of various phases of thermodynamic systems in terms of trajectories and RG flows in Hamiltonian space. In addition, 
it provided a scheme, to calculate, in principle, the behavior in the vicinity of critical and fixed points, by linearizing 
the flow equations. The method has proved to be powerful in providing an understanding of uniform (translationally 
invariant) systems, an exact implementation of the Wilson method in systems with quenched random disorder has been 
relatively restricted. Thus, for example, while the method has been successfully applied to the Anderson localization 
problem in (2 + e) dimensions,EI generalization to include electron-eleciron interactions^] have had limited success, 
i_i ' because of the complications such as the existence of local momentsu'El which are left out in such long- wavelength 
approaches to the problem. Another class of problems where a large amount of work has begn doneEl is based on 
the random antiferromagnetic (AF) spin- 1/2 chain in one dimension, where it has been showro that a perturbative 
real space renormalization group (RSRG) scheme becomes asymptotically exact, as has been explicitly -yerificd in a 
■ number of cases, such as the random XY chain, where it can be mapped on to a model of free fermions.S 
| Another issue where less attention has been given, is the range of validity of the asymptotic results of the linearized 
£T} i RG equations. In most models of uniform systems, a few iterations of the RG equations leads to the vicinity of the 
appropriate fixed point. Consequently, a study of the RG trajectory from the region in the vicinity of an critical 
fixed point to a stable fixed point, usually yields a relatively-, quick crossover behavior. However, in a study of a 
three dimensional model related to the ID random AF chainja it was found that the system did not settle down to 
a fixed-point-like behavior over several orders of magnitude in temperature^-4nstead there was a slow, logarithmic in 
temperature, evolution in the behavior. Indeed, it has been claimed recently£L9 that the true low temperature behavior 
of the model is likely different. Unfortunately, such a change cannot be experimentally investigated, because smaller 
. terms in the full Hamiltonian (such as hyperfine coupling to nuclear spins), and dipolar couplings, which have been 
neglected in the idealized model, will modify the behavior at such (exponentially low) temperatures. Nevertheless, 
it raises the issue that in models with quenched disorder, the true asymptotic fixed point behavior, may not be 
O ■ practically reached, at least in some cases, making "relevant" (in the RG sense) irrelevant in any practical sense. It 
is therefore of considerable interest to examine models with quenched disorder where crossovers between short length 
scale and long length scale behavior can be studied in some detail. 

One model, which ha^ recei | ye | d t ptpngiy c attention in recent years, is the random antiferromagnetic-ferromagnetic 
(AFM-FM) spin chain OH^IlSEallSOH The study of such systems has been motivated by the discovery of novel 
one-dimensional (ID) random spin systems. One example of such materials is Sr3CuPti_ x Ir x 06,c3 alloy of the 
pure compounds SraCuPtOg (antiferromagnet) and SraCuIrOg (ferromagnet). Such a system is interesting due to 
the interplay between quantum fluctuations and disorder; the latter can not be treated as a perturbation, as it is 
often found to change the low-energy spectra dramatically or even to destabilize the pure phases. Theoretical and 
numerical methods used to study random antiferromagnetic-ferromagnetic spin chaias have been largely inherited from 
the study of random antiferromagnetic (AEM) chains (with AFM coupling only).LrE3Eil In particular, the real-space 
renormalization group (RSRG) approachEjuj revealed a Curie-like temperature dependence in the susceptibility of a 
random AFM-FM spin- 1/2 chain at low temperatures. In the RSRG picture, the Curie behavior occurs because spins 
correlate, due to the existence of ferromagnetic couplings, to form clusters whose average size and effective spin grow 
in a random walk fashion at low temperatures. Therefore, the random AFM-FM spin- 1/2 chain belongs to a different 
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universality class from that of the random AFM spin-f/2 chain, whose ground state is known as a random singlet 
phase, in which singlets can be formed over large distance. The RSRG r esult s h ave been supported by various other 
numerical techniques, such-as quantum Monte Carlo (QMC) simulationE3E-0t3 and density matrix renormalization 
group (DMRG) methodMN 

In the RSRG method, one decimates the spins that are coupled by the strongest coupling in the system, and 
renormalizes the couplings among the remaining spins perturbatively. In this procedure new effective couplings and 
effective spins may also be generated. This procedure, which in general lowers the overall energy scale of the problem 
at hand, is repeated so that excitations at lower and lower energy scales are probed. Its quantitative accuracy 
relies on the presence of strong randomness, as in the case of strong randomness the strength of the strongest and 
typical couplings are well separated, ensuring the accuracy of a perturbative calculation. In real systems, however, 
the randomness can be very weak, especially when it is not introduced intentionally. In this case the RSRG is not 
expected to be quantitatively reliable in its early stages. It has been arguedBli-3 that the strength of disorder grows 
as the energy scale is lowered, thus in the low-energy (or equivalently, low-temperature) limit, the RSRG results 
become correct or even asymptotically exact. Nevertheless in this case one expects the thermodynamic properties of 
the system to be dominated by the physics of the pure chains at high temperature, and crossover to the randomness 
dominated regime (where the RSRG results apply) at very low temperatures. 

In the present work we study the thermodynamic properties of the random AFM-FM chain, with mostly uniform 
FM coupling and a small concentration of (impurity) AFM bonds, using the modified spin-wave method. We also 
compare the results of modified spin- wave method with those of exact diagonalization in small systems to demonstrate 
the accuracy of the former. Our motivation comes from the following considerations. First of all, as discussed above, 
when the concentration of the impurity bonds is small, the randomness is rather weak and the RSRG results are 
not reliable at high temperatures. The modified spin-wave method allows us to treat the high and low temperature 
regimes on equal footing, and in particular, address the crossover between these two regimes. This makes it possible 
to compare theory with experiments for the entire temperature range, when the randomness is weak. Secondly, the 
modified spin-wave method has so far been applied to study pure spin models only. While very successful in those 
cases, it has not yet been used to study random spin systems. By working out a number of technical issues that 
one faces when applying it to random spin problems, and demonstrate its accuracy in the present problem, we lay 
the ground for the application of this powerful method to other random spin problems. Thirdly, as discussed in the 
opening paragraphs, it is of general interest to study models where the short distance (high temperature) and long 
distance (low temperature) behavior are controlled by different fixed points in the renormalization group sense, and 
address the crossover between the two limiting cases quantitatively; our model is such an example, and we hope our 
study will stimulate future research on this important issue. 

Our results are summarized as follows. We find that the spin susceptibility of our system follows that of a pure 
ferromagnetic chain at high temperatures, and crosses over to 1/T dependence at low temperatures predicted by 
the RSRG, with a coefficient that agrees with the RSRG result essentially exactly. We also find that the crossover 
temperature depends mainly on the impurity bond concentration, while the width of the crossover depends mainly on 
the impurity bond strength; we determine their dependences semiquantitatively. We also demonstrate that the results 
of the modified spin-wave method is asymptotically exact in certain limiting cases, and quite accurate in the entire 
temperature range by comparing them with exact diagonalization results in the systems. 

The paper is organized as follow s. I n Sec. II, we describe the model we study, as well as the qualitative physics 



expected from the model. In Sec. Ill, we introduce the modified spin- wave theory for the random AFM-FM spin 



chains. We present the numerical results of the theory in Sec. IV, with emphasis on the effects impurity bond strength, 
impurity concentration, as well as finite-size effects of our numerical study on the spin susceptibility. Finally, our 
results are summarized in Sec. Readers interested in the comparison of the modified spin-wave theory to exact 
diagonalization results will find details in Appendix |A] and [b] for simple cases like single ferromagnetic segment and 
two coupled segments. 

II. THE QUALITATIVE PICTURE 

We consider the following Hamiltonian for a random spin- 5* chain 

H = Y,JiSi-Si+U (1) 

i 

where the coupling strength Ji is randomly chosen according to the distribution 

P(J) = (I - p)S(J + J F ) + p5(J - J AF ). (2) 
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Jf and Jaf are the coupling strengths of the ferromagnetic and antiferromagnetic bonds, respectively, p is the 
concentration of the antiferromagnetic impurity bonds. We are primarily interested in the dilute doping regime, i.e., 
p <C 1, so that we have long ferromagnetic spin segments, separated by antiferromagnetic bonds. The number of spins 
in a spin segment (N s ) is random, with an exponential distribution 

P(N s )=p(l-p) N °~ 1 . (3) 

Therefore, the average number of spins per segment is (N s ) = 1/p. From real space renormalization group studiesSE], 
we understand that a strong antiferromagnetic bond tends to lock the two adjacent spins into a singlet (which is 
inert at temperatures of interest), and the singlet thereby generates a weaker effective AF coupling between the two 
spins next to the singlet (see Appendix ^). Therefore, for simplicity, we consider weak antiferromagnetic coupling 
(Jaf < Jf) on ly without losing generality. 

The temperature dependence of magnetic susceptibility x per spin is qualitatively known. Wc assume the Zccman 
term H' = —^H z ^2,Sf and set /j, = ks = 1 for simplicity. At high temperatures (T 3> Jf), the spins behave 
independently. Therefore, we expect an ordinary Curie contribution of 

XHT = 3T , (4) 

from each spin. When T < Jp, the spins within each ferromagnetic segment become correlated. TakahashS obtained 
the following low temperature expansion for spin- 1/2 ferromagnetic chain (of infinite length) 
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using the Bethe-ansatz integral equations. The factor of 1/4 is added here due to an extra factor of 2 in the Zeeman 
energy terms in the original paper. Therefore, susceptibility of each segment rises as temperature decreases, crossing 
over from 1/T to 1/T 2 . In the case of Jaf = 0, X crosses over back to Curie behavior at lower temperatures, when 
each segment of finite length acts like a block spin with frozen internal excitations. When T drops below the typical 
spin-wave excitation gap of a typical segment 

A sw ~J F s(^j\ (6) 

we expect 

cbs _ 1 N„S(N.S + 1 



1_ v N S S(N S S + 1) 



XBS t 3T 

s 

Here N s is the typical length of a ferromagnetic segment. In the thermodynamic limit, we obtain the Curie constant 
for locked block spins to be 

S 2 2-p S 

CBS = T— + r (8) 

using segment-size distribution in Eq. |^. At p — 0.1, we obtain cbs — 7/4 for S = 1/2. cbs decreases with increasing 
impurity concentration p and is always (unless p = 1, i.e. in the pure antiferromagnetic case) greater than the high 
temperature Curie constant (in Eq. ^) 

S(S + l) 1 1 
cht = - -, for S=~, (9) 

expected at high temperatures. 

When the end spins of neighboring segments are coupled with weak antiferromagnetic coupling Jaf, the block spins 
(originating from the ferromagnetic segments) interact with weak effective couplings. For two segments with N\ and 
N2 spins, we obtain the effective coupling 
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by projecting the coupling to the block-spin snace using Wigner-Eckart theorem. Not surprisingly, long segments 
are weakly coupled. In the RSRG approach,E3 one starts with two strongly coupled segments and decimates the 
segments to a S' — \N\ — Ns\S spin, and renormalize its couplings with neighboring segments. The spin- wave 
excitations within each segment, considered in our modified spin- wave approach, are in general neglected in the RSRG 
approach, which is, therefore, valid only at low enough temperatures. By including these spin-wave excitations, we 
can not only demonstrate that the Curie susceptibility (predicted by RSRG) in a AFM-FM spin-1/2 chain at low 
temperatures is indeed the result of correlated ferromagnetic segments, but estimate the temperature scale at which the 
low temperature Curie behavior occurs as well. Since the renormalized couplings may become ferromagnetic, clusters 
start to grow in a random walk fashion as segments couple to neighboring segments as we lower the temperature. As 
in a typical random walk problem, we expect the average spin Si and the average size I of such clusters of spins satisfy 

Si ~ r 1/2 , (n) 

so that as far as magnetic susceptibility per spin \ is concerned, the scaling behavior of Sf and I cancels, leading to 

XiT= ^ = lS.I. (12) 
T 3T I T y ' 



By statistical analysisjll one finds the low-temperature Curie constant to be 

clt = -7, . (13) 

3 p 

in the large cluster-size limit. This result coincides with the classical result for the susceptibility per spin at low 
temperaturesOli3 Compared with Eq. ^, we find that antiferromagnetic impurity bonds reduce the low-temperature 
Curie constant by 

S 2 S 

Ac = c B s-c LT = — + -. (14) 
3p 3 

On the other hand, for small p (or low concentration of AF bonds), we have clt > cht] this simply reflects the fact 
that the dominant interaction in the system is ferromagnetic, which enhances spin susceptibility. Thus in the cases 
we are interested in, we have 

cht < c LT < cbs- (15) 

We should note here that unless Jaf is very small, there may not be a temperature range in which \ ~ Xbs = 
cbs/T] this is because the neighboring segments can start developing correlations due to Jaf before the intrasegment 
spin-wave excitations are frozen out by temperature. We will actually study mostly this case in this work, and 
focus mainly on the crossover from the pure ferromagnetic chain behavior directly to the asymptotic low temperature 
behavior with clt as the Curie coefficient. We do this partly because in real systems Jaf is often comparable to Jf, 
and partly for the sake of simplicity. Our method, however, is capable of handling very small Jaf and the regime 
with cbs being the Curie coefficient. 



III. THE MODIFIED SPIN- WAVE APPROACH 



The conventional spin- wave theory gives the exact spin- wave spectra of Heisenberg ferromagnets at T — 0. At finite 
temperature however, the theory leads to difficulty in one-dimensional quantum ferromagnets, the number of spin 
waves diverges when external magnetic field goes to zero. This difficulty has its root in the Mermin- Wagner theorem, 
which dictates that there cannot be long-range order in ID at finite temperature. Bv imposing a constraint that the 
total magnetization be zero (thus fulfilling the Mermin- Wagner theorem) , TakahashO introduced a chemical potential 
(equivalence to a uniform magnetic field) for the spin waves and successfully obtained the low-temperature properties 
of the ferromagnets in one dimension, as well as in two dimensions. The beauty of this modified spin-wave theory is 
that it gives correct asymptotic thermodynamic behavior at both large and zero temperature limit, as well as fairly 
accurate crossover behavior at intcr-eaecLiate temperatures. The modified spin- wave theory have further succeeded 
in two-dimensional antiferromagnetscZIES and one-dimensional ferrimagnets.ESLFor one-dimensional antiferromagnets, 
the modified spin-wave theory predicted a gap in the spin-wave spectrum,E3 consistent with the Haldane gap in 
integer-spin chains.EJ 
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In this paper, we extend its applicability to ferromagnetic spin chains with dilute antiferromagnetic impurities. In 
the dilute impurity limit, a spin chain consists of ferromagnetic spin segments, coupled antiferromagnetically. The 
application of the modified spin- wave theory to the random antiferromagnetic-ferromagnetic spin chains is, therefore, 
mainly based on our observation that ferromagnetic spin segments can be well described by the modified spin-wave 
theory (see Appendix |A|) . We found that, for magnetic susceptibility of a finite spin chain, the modified spin- wave 
theory gives noticeable difference from the exact diagonalization result only at two temperature crossovers. The 
crossover at higher temperatures describes the switching on of the couplings between individual spins, which is of 
little interest in the competition between different segments. On the other hand, the maximum error of the modified 
spin-wave theory is roughly 5%, for a chain of up to 14 spins, near the lower-temperature crossover when the spin 
chain becomes locked into one block spin. For a finite spin chain, the translational symmetry is broken. Therefore, 
one, in principle, needs to introduce a local chemical potential (or a local field) for each spin to ensure the magnetic 
moment of the spin be zero. Fortunately, we can apply a periodic boundary condition to each ferromagnetic spin 
segment, so that only one chemical potential is needed throughout the segment. We emphasize that the use of periodic 
boundary condition reduces the computational complexity without changing the essential physics. This is particularly 
true when we have very dilute impurities so that the average ferromagnetic segment length is large. 

For the nearest-neighbor Heisenberg Hamiltonian, Eq. ^, we apply Holstein-Primakoff transformation 




= J2S-afaial, 
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for segments with odd index s, and 
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for segments with even index s. This way, we introduce a distinct species of boson for each segment. The Hamiltonian, 
in the linear spin-wave approximation, becomes 



H = E Q + J F Sj2J2{ a ? a t + a &i a °+i-^<+i- a &i a i 



s=l i=l 

M s -I 



JafS ^ { a N a -l a N s -l + a l 



l N s -l u 'l 



l N s -l I i 



(18) 



where M s is the number of ferromagnetic segments, and N s the number of spins in the s-th segment. The total 
number of spins is, therefore, N = 53 s =i Note (N s + 1) = 0, reflecting the periodic boundary condition imposed 
on each segment. 

In the linear spin- wave approach, the quadratic Hamiltonian is soluble by a generalized Bogoliubov transformation 
for bosons, a generalization of Ref. ^J. The procedure can be written compactly in a matrix format. Denote the 
original boson operators by a vector 



_ / 1 1 1 2t 2f 2f \1 

x — (a 1 , a 2 , ■ ■ ■ , cl Ni , a-L , a 2 , • ■ • , a N2 , ■ ■ ■) 



The Hamiltonian can be written as 



H = const. + x'Hx, 



(19) 
(20) 
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where 



/ 2 
-1 

-1 



Tt = JpS 



-1 

-1 2 + Jaf/Jf Jaf/Jf 

Jaf/Jf 2 + Jaf/Jf — 1 

-1 2 



-1 



-1 

-12 + J afI Jf 



Jaf/ Jf 



Jaf/Jf 2 + Jaf/Jf 



V 



(21) 



/ 



is a N x iV matrix. The Hamiltonian matrix is basically a 3-diagonal matrix, with extra (-l)'s at the corners of 
each block (of a spin segment). The rest of the matrix elements (left unspecified in Eq. pl|) are zeros. Introduce the 
generalized Bogoliubov transformation: 



x = Vj, 

where V is a N x N matrix and 7 a vector of the set of boson operators (a's) 

/ll 1 2t 2t 2f \T 

7 = («i,a 2 , • • >«2>' ' ' > a W 2 > • ' ') > 

that diagonalizes Tt, i.e., 

V T HV = E, 

in the compact matrix format, where E is a diagonal matrix, 

^JVxiV — diag{e 1 , e 2 , • • • , e^, £ x ,e 2 , • • ■ , e^, • • ■ , e Wm j 



(22) 
(23) 
(24) 
(25) 



with its diagonal elements being the energies of the corresponding bosons after the Bogoliubov transformation. The 
boson commutation relation requires 

VPV T = P, (26) 
where P is a generalized para unit matrix (inheriting the notation from Ref. |2~i| ) of the form 

Ni N 2 

P NxN =diag{lX^~~ 1,-1, -V ••,-!,•■■) (27) 

For fermions, on the other hand, the fermionic commutation relations lead to the orthonormalization condition, which 
can assume the same equation as Eq. |2^, except that P becomes an N x N unit matrix. We need to find out a 
solution V for bosons, which simultaneously satisfies Eq. U and ^6|; this can be done with the help of generalized 
matrix diagonalization. Unlike fermionic diagonalization, the existence of such a solution is not guaranteed, unless 
spin segments are decoupled, i.e. Jaf — 0. Such unfortunate situations are, in general, associated with modes of 
zero energy, which then leads ta divergent number of bosons (or divergent magnetization). In the case of a single 
ferromagnetic chain, TakahashiES introduced a chemical potential (or a magnetic field) to overcome the difficulty. In 
this case, we need to introduce to the Hamiltonian a set of M s local chemical potentials, one for each segment, to find a 
solution. These chemical potentials are chosen such that the magnetization of each spin segment (equivalently, of each 
spin with periodic boundary condition) is zero. These constraints, which have the physical significance of restoring 
the rotational symmetry, ensure that the number of bosons be finite at finite temperatures. The thermodynamic 
quantities can be calculated from the excitation energies e's and the transformation matrix V. In particular, magnetic 
susceptibility per site x can be written as 



1 



X 



3TN 



s=l t=l 



where hf is the occupation number for the boson correspond to operator af and excitation energy sf 

1 



l/T 



(28) 



(29) 
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IV. RESULTS FOR MAGNETIC SUSCEPTIBILITY 



Figure |l| shows the static magnetic susceptibility \ per spin of a 600-spin ferromagnetic chain, with 59 antiferro- 
magnetic impurity bonds, which divide the chain into 60 ferromagnetic segments, i.e., p = 0.1. The impurity bonds 
(Jaf = 0-5) are put randomly with the constraint that the total magnetization of the corresponding classical ground 
state is zero. This guarantees that \ drops to zero at zero temperature. The susceptibility curve shows three temper- 
ature regimes. For T/Jp > 3, result can be fit to a high-temperature Curie law (Eq. ||), suggesting that spins behave 
independently. For 0.2 < T/Jp < 3, \ rises with decreasing temperature, following Eq. |[ This implies that spins 
start to correlate, forming independent ferromagnetic segments. The antiferromagnetic coupling is still too weak to 
affect thermodynamics of the random spin chain at these temperatures. Below T — 0.2Jp, x is m good agreement 
with 

X = — , CLT =4' ( 30 ) 

the low-temperature Curie behavior expected by the RSRG approach (Eq. [l^ for p = 0.1). We point out that 
the Curie constant of independent block spins (Eq. |J) is expected to be cbs — 7/4 for p — 0.1. Therefore, the 
low-temperature Curie constant is indeed reduced by antiferromagnetic couplings between ferromagnetic segments. 
We note that x deviates from Eq. ^ below T = 0.02 Jp. We believe this is an artifact due to finite size (we have 
only 600 spins), since there are no more segments to decimate, x bends down as we specially require that x S° to 
zero at zero temperature. Since the number of iterations required to find a solution for the generalized Bogoliubov 
transformation grows significantly at low temperatures (nearly 1000 iterations at T — 0.005), we were unable to 
approach temperatures much lower than T — 0.01 or to calculate many samples on our regular workstations. It 
is worth pointing out that we can identify the crossover temperature from the ferromagnetic spin-chain physics to 
random spin-chain physics as the cross point of Eq. ^0| and |^. In the dilute doping limit (small p), we have 



T x = Jf: 



0.50 + 3.04 J +0 



1 — p \1 — p 



(31) 



For p = 0.1, T x =0.168J F . 

Figure |^ shows the susceptibility per spin x °f spin chains of 20, 30, and 40 ferromagnetic segments (i.e., with 
19, 29, and 39 antiferromagnetic impurity bonds, respectively). The average ferromagnetic segment length is fixed 
at 8 spins (p — 1/8), so that the lengths of the corresponding spin chains are 160, 240, and 320, respectively. The 
antiferromagnetic coupling Jaf is chosen again to be 0.5. x is averaged over 10-20 random realizations depending 
on size. For 0.2 < T/Jp < 1, x roughly follows the low-temperature (compared with Jp) expansion of a single 
ferromagnetic chain (Eq. ^J). From now on, we neglect the thermodynamics at T > Jp (although it is capable of being 
calculated within our theory), since it is trivial and not the interest of this paper. Below the crossover temperature 
T x ~ 0.2, x deviates from the pure chain susceptibility, and obeys a Curie law. The range of the Curie susceptibility 
extends to lower and lower temperature when we have more segments (from 20 to 40 in Fig. ^), therefore longer chain, 
for a fixed impurity concentration. At N s = 40, this range extends beyond one decade for p = 1/8. We point out 
that the Curie constant obtained from Fig. || (c = 0.67) is roughly 15% larger than predicted by Eq. [l3|. This should 
be viewed as a mixture of both the RSRG- predicted behavior and the ferromagnetic-chain physics. The reason will 
become clear later after we explore the impurity concentration dependence of the Curie behavior. 

The antiferromagnetic coupling strength Jaf determines the temperature scale at which the ferromagnetic segments 
couple to their nearest neighbors. The smaller Jaf is, the lower the temperature at which RSRG results are valid 
becomes, leading therefore to an increasingly wider crossover from the physics of a pure ferromagnetic chain. Figure |3| 
shows the sample averaged magnetic susceptibility per spin x for Jaf — 0.3 and 0.5 for spin chains of 30 ferromagnetic 
segments. The average number of spins in each segment is 12 spins (p = 1/12). The dashed line corresponds to a 
Curie constant c = 11/12, expected by Eq. [l3] for p = 1/12. Although the system size (360 spins) is probably not 
large enough, one can nevertheless identify the Curie behavior for Jaf = 0.5 around T — 0.05 and for Jaf = 0.3 
around T = 0.02. Note that x f° r Jaf = 0.3 is greater than that for Jaf — 0.5 at low temperatures. This is 
expected because decimated spin segments, in the RSRG language, contribute less to the magnetic susceptibility 
than the sum of original segments, and because the decimation happens at lower temperatures for smaller Jaf- For 
very small Jaf, we expect the system crosses over from a ferromagnetic regime at high temperatures, first, to an 
(almost) independent block-spin regime at lower temperatures, characterized by a Curie constant as in Eq. |^. There 
is another crossover, at very low temperatures (dependent on Jaf), from the independent block-spin regime to the 
RSRG regime, characterized by a smaller Curie constant as in Eq. [l3| 

We have also studied the dependence of the magnetic susceptibility on the impurity concentration, or equivalently, 
average number of spins in each segment. Qualitatively, the larger number of spins within each segment increases the 
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FIG. 1: Static magnetic susceptibility per spin (x) of a single random Heisenberg spin chain of 600 spins. The spins are coupled 
with their nearest neighbors by ferromagnetic couplings (Jf), with 59 antiferromagnetic impurity couplings (Jaf = 0.5Jf), 
so that the spin chain can be viewed as 60 ferromagnetic segments of average length 10 (spins). On the log- log scale, the two 
dashed lines with slope unity are exact results of the ordinary Curie susceptibility of independent spins at high temperatures 
(Eq. ^|) and the low-temperature Curie susceptibility expected by the RSRG approach (Eq. |l3| ). The dotted line is the low 
temperature expansion of the magnetic susceptibility (Eq. ^) of a pure spin- 1/2 ferromagnetic chain (of infinite length), which 
is proportional to 1/T 2 in the low temperature limit. 
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FIG. 2: Sample averaged magnetic susceptibility per spin for spin chain of 160, 240, 320 spins, which consists of 20, 30, 40 
ferromagnetic segments, respectively, coupled by impurity antiferromagnetic bonds with Jaf = 0.5 Jf- The dotted line is the 
low temperature expansion of the magnetic susceptibility (Eq. ||) of a spin-1/2 ferromagnetic chain (of infinite length). The 
dashed line is a fit to the apparent (see text for detail) low-temperature Curie susceptibility. 
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FIG. 3: Sample averaged magnetic susceptibility per spin for spin chains of 360 spins in 30 ferromagnetic segments, coupled 
by antiferromagnetic bonds with Jaf = 0-SJf and 0.5 Jf- The data is expected to agree with the low-temperature Curie 
susceptibility (the dashed line) in the RSRG method. The temperature range of the agreement is expected to be higher for 
stronger antiferromagnetic impurity bonds (Jaf = O.BJf). The dotted line is the low temperature expansion of the magnetic 
susceptibility (Eq. |^) of a spin- 1/2 ferromagnetic chain (of infinite length). 



magnetic susceptibility per spin, since \ ~ for a ferromagnetic A^-spin segment (for large N s ). As a consequence, 
the Curie constant is expected to increases with the average number of spins in each segment, or more precisely, 
as in Eq |l3|. This trend is shown in Fig. |[ where we plot the averaged susceptibility per spin for 30-segment spin 
chains, with 8, 12, and 16 spins per segment on average, respectively. The three dashed lines are the low-T Curie 
behavior the RSRG approach expected (Eq. |l3[ ). As is pointed out in earlier paragraph, our data is in significant 
disagreement with Eq. 1^ for p = 1/8. For p = 1/16, one finds that x follows Eq. [l^ from T = 0.015Jp to 0.04, 
the regime antiferromagnetic impurity bonds are dominating. The Curie regime shifts to 0.02 < T/Jp < 0.08 for 
p = 1/12. This is consistent with Eq. that longer segments have relatively weaker effective intersegment coupling. 
The trend suggests that for p — 1/8 the Curie regime would occur above T > O.lJp, which coincides with the 
regime where intersegment coupling is important. Therefore, we believe the discrepancy between our data and Eq. p"3] 
reveals a collective contribution from the ferromagnetism with each segment and the antiferromagnetism between 
segments. This demonstrates that inappropriate parameters can kill (or lead to the wrong) RSRG physics. We note 
that the crossover regime (deviating from both Eq. ^ and the 1/T-law) becomes wider for smaller p. This is again 
consistent with the translation between low concentration and weak effective intersegment coupling, which, we have 
demonstrated, leads to wider crossover regime. 



V. CONCLUSION 



In the study of the random AFM-FM (nearly ferromagnetic) spin-1/2 chain, the modified spin- wave theory al- 
lows us to obtain the thermodynamic properties over a broad temperature range (on logarithmic scale). We have 
studied systems as large as 600 spins, with up to 60 ferromagnetic segments, for temperature ranges of over three 
orders of magnitude. At very high temperatures, the spins are independent; correlation starts to develop within 
each ferromagnetic segments as temperature drops, giving rise to characteristics of a pure ferromagnetic chain. At 
low temperatures, the observation of the Curie-like behavior reveals the regime where ferromagnetic segments are 
correlated into clusters, whose size grows with decreasing temperature. Therefore, we can study, fairly systematically, 
in the modified spin- wave theory the crossover between the high-temperature Curie behavior (independent spins) and 
low-temperature Curie behavior (correlated spin segments). 
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FIG. 4: Sample averaged magnetic susceptibility per spin for spin chains of 240, 360, and 480 spins, each containing 30 
ferromagnetic segments, coupled by antiferromagnetic bonds with Jaf = 0.5 Jf- The dotted line is the low temperature 
expansion of the magnetic susceptibility (Eq. |^) of a spin-1/2 ferromagnetic chain (of infinite length). The dashed lines are the 
exact results of the low-temperature Curie susceptibility expected by the RSRG approach for the three impurity concentrations. 

System size is important in observing the low-temperature Curie behavior, since it determines whether the RSRG 
scaling regime can be reached. Our study has found that the modified spin-wave theory gives results in accordance 
with the RSRG argument for chains of several hundred spins (after averaging over random bond distribution). The 
study of the effects of impurity bond strength shows that the low-temperature Curie regime can be consistently pushed 
to lower temperatures by decreasing the antiferromagnetic impurity bond strength. 

In the dilute doping regime, we have studied the effects of impurity concentration on the low-temperature ther- 
modynamics. The results are again in fairly good agreement with the RSRG prediction. The exception happens 
at higher concentration (p = 1/8), where we believe the intrasegment ferromagnetic coupling and the intersegment 
antiferromagnetic coupling take actions at about the same temperature regime, leading to a mixed behavior which 
appears to look like a different Curie behavior. Further explorations using other methods would be very useful to 
confirm our explanation. 

In our study, both weak impurity bonds and low impurity concentration give rise to a wide crossover from the 
ferromagnetic spin-chain regime to the random spin-chain regime. Although the modified spin-wave theory limits 
us to only a portion of the whole random spin-chain space, we believe such wide crossover can be generic and it 
can prevent the real R-SftG behavior from being observed at reasonably low temperatures. In fact, the RSRG study 
of random spin chainsLj has pointed out the possibility of having a crossover region of more than five orders of 
magnitude before the tnie RSRG scaling regime can be reached. A quantum Monte Carlo study of the random 
AFM-FM spin-1/2 chair£j has confirmed that a region of weakly interacting spin segments can exist, whose crossover 
to the low-temperature scaling regime is marked by a (somewhat small) peak in the specific heat. 

The main difficulty of directly comparing our results with experiments seem to be that quasi-lD Heisenberg ferro- 
magnetic systems are rare since there usually exist interchain couplings which lead to three-dimensional (3D) ordering 
at low temperatures. For SraCuIrOg, long-range (presumably 3D) order develops at 20.1 K without signature of 
ID magnetism.Ej Recently, however, effects of magnetic impurities en the quasi-lD ferromagnetic spin chain have 
been reported for organic radical alloy, (p-CDpOV)i_ 2 ,.(p-BDpOV) 2; ,E3 whose interchain-intrachain coupling ratio is 
as small as J'/ J — 3.7 x 10~ 3 . This may open a door where the low-temperature Curie susceptibility of the RSRG 
origin can be observed in ferromagnetic bond dominated random spin chains. 
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APPENDIX A: MODIFIED SPIN- WAVE THEORY FOR FERROMAGNETIC SPIN SEGMENTS 

In Takahashi's original paper@, the modified spin-wave theory is aimed at ferromagnetic spin chains with infinite 
length. Nevertheless, the theory can be easily applied to a finite-size ferromagnetic spin chain with periodic boundary 
conditions. This only changes the continuous spin- wave spectrum to a discrete spectrum. Therefore, for an iV-spin 
ferromagnetic chain, we can write the self-consistent equations for the single (global) chemical potential ijl: 



NS = ^n fe , 



"* ~ e [eW-p]/T_!' ( A1 ) 



where the discrete spin- wave energy spectrum is 



2tti 

e(k) = 2 J SCI -cosfca), k= , i = 0, 1, N - 1. (A2) 

Na 



The magnetic susceptibility per spin can be expressed as 



^5> fc (n fe + l). (A3) 



x -- 3TN 

k 

Analytical results, which turn out to be exact, can be obtaine d fo r both high and low temperature limit. In the 
high temperature limit, = S is a constant. Self-consistent Eqs. |Al] give 

IM = Tin (l + S- 1 ) . (A4) 

Therefore, we have 

(A5) 

which is the exact result for inde pend ent spins. In the low temperature limit, only k — mode can be excited, i.e. 
hk = NSSk,o- Self-consistent Eqs. Al give 



/x = Tin (l + N~ 1 S~ 1 ) . (A6) 

Therefore, we have 

x = ^>, (A7) 

which is the exact result for the ferromagnetic ground state with maximum spin. 

Figure [3] shows x f° r ferromagnetic spin- 1/2 chains of finite length N — 4 and 14. The modified spin- wave theory 
results agree very well with exact results for both sizes. The largest descrepency (less than 15%) occurs around 
T = Jp, at which independent spins start to correlate. This regime is of less interest for the purpose of this paper. 
In both the high temperature limit (T > 10 Jp) and the low temperature limit (T < Asw), the modified spin- wave 
theory gives exact results. This is a very impressive result since the spin-wave theory is in principle developed around 
T = 0, describing the low temperature excitations of the ferromagnetic spin chain. In addition, we are able to show 
that boundary conditions have only weak effects on the rapidly changing susceptibility. In particular, behavior in the 
high temperature limit and the low temperature limit remains unchanged. 
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FIG. 5: Exact diagonalization result (solid line) and the modified spin-wave theory result (dashed line) of the magnetic 
susceptibility per spin for (a) N = 4 and (b) N = 14 ferromagnetic spin segments (S = 1/2) with periodic boundary condi- 
tions. The two dotted lines are the high-temperature Curie susceptibility of independent spins and the low-temperature Curie 
susceptibility of the spins when frozen into a block spin. 



APPENDIX B: COUPLING BETWEEN TWO FERROMAGNETIC SPIN SEGMENTS 



The problem of two ferromagnetic spin segment coupled by an antiferromagnetic bond is a special case of more 
general random spin chains. For small systems, we can compare our modified spin- wave theory with exact diago- 
nalization. Figure ^ shows the susceptibility for two ferromagnetic segments (Jp = 1) coupled with a Jaf = 0.01 
antiferromagnetic bond. The numbers of spins in the two segments are Ni = 6 and N 2 = 7. Three regimes where the 
susceptibility obeys 1/T-law are well-defined in Fig. ^. The corresponding pictures from high to low temperature are 
independent spins, two independent spin segments, completely locked spins, respectively. The susceptibility can be fit 
very well to the Curie law obtained from the simple pictures. In these Curie regimes, the modified spin-wave theory 
gives exact results. Discrepancies between the modified spin-wave results and exact results can only been noticed 
at the two crossover regimes between them. Crossover temperature scales can be roughly estimated by the smaller 
spin-wave gap Asw of the two segments and the effective antiferromagnetic coupling J' between the two segments, 
as marked by arrows in Fig. 

In the RSRG scheme, one introduces an effective antiferromagnetic coupling J' between the two segments, which 
replaces the coupling between the two neighboring end spins. By inspecting the energy spectrum of the two-segment 
system in Fig. ffl, one can explore the validity of the replacement. At low enough temperatures (T < Asw), the 
finite spin-wave excitations are frozen within each segment. Therefore, only k = modes are relevant when the two 
segments couple with the effective J'. This results a group of lowest energy levels (S = 1/2, 3/2, 13/2), which are 
separated from higher energy levels by a energy gap (~ Asw)- The gap is expected to become comparable with the 
width of the k = energy levels, when Jaf is close to Jf, when RSRG argument becomes less accurate. For small 
Jaf, these k = modes obtained from the exact diagonalization are in good agreement with the effective levels in the 
RG picture (the inset of Fig. ^). The difference in energy levels given by the two methods grows with increasing Jaf- 
Nevertheless, the agreement in the susceptibility obtained from the two methods still persists even when Jaf ~ Jf, 
as long as T <C Asw, which poses a temperature restriction on the RSRG scheme. 



1 K. G. Wilson and J. Kogut, Phys. Rep. 12, 75 (1974). 

2 F. Wegner, Z. Phys. B 25, 327 (1976); Z. Phys. B 35, 207 (1979). 



13 



100000 
10000 
1000 
100 
10 

1 

0.1 
0.01 



0.001 



N 1 =6, N 2 = 7, S = 1/2 



J F = 1, J AF = 0.01 



Exact 
Spin-wave 



S(S+1) 
3T(N 1+ N 2 ) 



N 1 S(N 1 S+1)+N 2 S(N 2 S+1) 
3T(N 1+ N 2 ) 



J' 




1e-06 1e-05 0.0001 0.001 0.01 

T/J c 



0.1 



10 



100 



FIG. 6: Exact diagonalization result (solid line) and the modified spin-wave theory result (dashed line) of the magnetic sus- 
ceptibility per spin of two spin-1/2 ferromagnetic segments (JVi = 6 and A^2 = 7) coupled by a Jaf = O.OUf antiferromagnetic 
bond. The three dashed lines describe the Curie behavior in the following three regimes. (1) At high temperatures, spins are 
independent. (2) At intermediate temperatures, spins within each ferromagnetic segments are locked into a block spin with 
no internal excitations. The two block spins interacts with each other through the effective antiferromagnetic coupling (J') 
generated by the antiferromagnetic bond (Jaf)- (3) At low temperatures, the two block spins are locked into one S = 1/2 spin. 
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FIG. 7: Energy spectrum of the two ferromagnetic segments described in Fig. ^| (JVi = 6 and N2 = 7 spin segments coupled 
by a Jaf = 0.01 antiferromagnetic bond). The lowest energy manifold (k — modes) of the exact diagonalization result (solid 
lines), shown on an amplified scale in the inset, is compared with the energy spectrum of the two block spins coupled by the 
effective coupling J' (Eq. ^o| ). 
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